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We present a theoretical study of the density of states and supercurrent in diffusive 
superconductor-normal metal-superconductor (SNS) junctions. In particular, we study the influ- 
ence on these two equilibrium properties of both an arbitrary transparency of the SN interfaces and 
the presence of spin-flip scattering in the normal wire. We show that the minigap that is present in 
the spectrum of the diffusive wire is very sensitive to the interface transmission. More importantly, 
we show that at arbitrary transparency the minigap replaces the Thouless energy as the relevant 
energy scale for the proximity effect, determining for instance the temperature dependence of the 
critical current. We also study in detail how the critical current is suppressed by the effect of spin-flip 
scattering, which can be due to either magnetic impurities or, under certain circumstances, to an 
external magnetic field. Our analysis based on the quasiclassical theory of diffusive superconductors 
can be very valuable to establish quantitative comparisons between experiment and theory. 



I. INTRODUCTION 

When a normal metal (N) and a superconductor (S) 
are brought together, their mutual interaction results in 
the modification of their electronic and transport proper- 
ties. In particular, the normal metal may acquire genuine 
superconducting properties such a gap in the density of 
states or the ability to sustain a supercurrent. This ef- 
fect, known as proximity effect, was first discussed by de 
Gennes 1 - in the 1960's and in recent years it has been 
extensively studied in diffusive hybrid nanostructures^ 
Many equilibrium 3 -!! and transport properties 5 -^ of diffu- 
sive SN systems are now well understood, which is par- 
tially due to the impressive predictive power of the qua- 
siclassical theory of superconductivity for diffusive sys- 
tems, which is summarized in the Usadel equations Z 

The proximity effect is mediated by Andreev reflec- 
tions!* In this tunneling process an electron coming from 
N with energy e below the superconducting gap A is con- 
verted into a reflected hole, thus transferring a Cooper 
pair to the S electrode. The time-reversed states in- 
volved in this process are coherent over a distance Lq = 
min(y hD/e, L$), where D is the diffusion constant of 
N and is the phase coherence length. This coher- 
ence may also be altered by interactions that break the 
time-reversal symmetry such as those induced by para- 
magnetic impurities or an external magnetic field. 

In this work we present a theoretical study of the den- 
sity of states and the supercurent in diffusive SNS junc- 
tions. These quantities nicely reflect the proximity ef- 
fect under equilibrium conditions. It was first shown 
by McMillan^ that a diffusive normal metal in contact 
with a superconductor can develop a gap in its elec- 
tronic spectrum, which is usually referred to as mini- 
gap. More recently the minigap has been studied by nu- 
merous authors, usually within the framework of the Us- 



adel equations . 10 ' 11 ' 12 ' 13 ' 14 From the experimental point 
of view, the appearance of a minigap has been tested in- 
directly in several tunneling experiments (see for instance 
Refs. \m flfl ] an references therein). 

On the other hand, the fact that a SNS junction can 
sustain a supercurrent is known since the first experi- 
ments performed with Pb-Cu-Pb sandwiches . 1 7 i 18 It was 
soon realized that the existence of a dissipationless cur- 
rent in these structures is possible due to the proxim- 
ity effect. 1 Later on, a more systematic experimental 
study of the critical current in these hybrid structures 
was carried out with the help of diffusive SNS micro- 
bridgesj 19 ' 20 The results of these experiments were de- 
scribed by Likharev,— who made use of the Usadel equa- 
tions in the high temperature regime (A <C &bT). A 
more general study of the Josephson effect in diffusive 
SNS junctions was made in Ref. [22]. More recently, Du- 
bos et aL— demonstrated that the full temperature depen- 
dence of the critical current of diffusive Nb-Cu-Nb junc- 
tions with highly-transparent interfaces could be quanti- 
tatively described by the quasiclassical theory. 

Most of the theoretical work done on proximity ef- 
fect has been focused in the case of either ideal (per- 
fectly transmissive) SN interfaces or in the tunneling 
limit , 23 ' 24 ! 25 with some notable exceptions. 2 ^ One of the 
two main goals of this paper is to study how the local 
density of states (DOS) and the supercurrent in diffusive 
SNS junctions is influenced by arbitrary transmission of 
the interfaces. This is an important issue, in particular, 
in order to be able to establish quantitative comparisons 
between theory an experiment, since in reality the mis- 
macht of material parameters leads to a broad range of 
transmission through the SN interfaces. In particular, 
we shall discuss the following issues: (i) how the trans- 
mission determines the magnitude of the minigap and, in 
general, the shape of the DOS in the normal wire, both 
in the absence and in the presence of a supercurrent and 
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(ii) how a finite transmission modifies the current-phase 
relation and the critical current of these junctions. Our 
results, based on the quasiclassical theory, show that the 
minigap, which is reduced as the interface resistance in- 
creases, is the energy scale that controls, in particular, 
the magnitude and temperature dependence of the crit- 
ical current. For ideal interfaces, this role is played by 
the Thouless energy ex = HD / L 2 , where L is the length 
of the normal metal. 

The second goal of our work is to study the role of spin- 
flip scattering in the local DOS and supercurrent of diffu- 
sive SNS junctions. This type of scattering, which can be 
induced by magnetic impurities or an external magnetic 
field, breaks the time reserval symmetry between the elec- 
trons in Cooper pairs and reduces the superconducting 
correlations ! 27 ' 28 Different author a 11 ' 14 ' 24 ' 29 ' 30 have stud- 
ied the effect of spin-flip scattering in the properties of 
SN structures. However, basic quantities like the super- 
current in SNS structures have, to our knowledge, not 
yet been addressed. We present in this work a detailed 
study of the dependence of the critical current of a dif- 
fusive SNS on the strength of the spin-flip scattering for 
arbitrary range of parameters (length of the normal wire, 
temperature and interface resistance). Our predictions 
can tested experimentally by measuring the critical cur- 
rent in the presence of magnetic field since, as long as 
the normal wire is a thin film, the field acts simply as a 
pair-breaking mechanism. 11 ' 28 ! 31 

The rest of the paper is organized as follows. In the 
next section we describe the general formalism, based on 
the quasiclassical theory for diffusive superconductors, 
that allows us to compute the DOS and the supercurrent 
in diffusive SNS junctions for arbitrary length, temper- 
ature and interface transparency. Section III is devoted 
to the analysis of the local DOS in the normal wire in 
different situations. In section IV we discuss the results 
for the dependence of the supercurrent on the interface 
resistance, temperature and strength of the spin-flip scat- 
tering. Finally, we summarize our main conclusions in 
section V. In the Appendix A we describe our analytical 
results for the critical current in the limit of weak proxim- 
ity effect and in Appendix B we include some numerical 
fits described in the previous sections. 



II. QUASICLASSICAL GREEN'S FUNCTIONS 
FORMALISM 



We consider the SNS junction represented schemati- 
cally in Fig. [TJ where N is a diffusive normal metal of 
length L coupled to two identical superconducting reser- 
voirs with gap A. We assume that the transport is phase- 
coherent, i.e. L <C L<f, and neglect the suppression of the 
pair potential in the S leads near the interfaces. Our 
main goal is to study how the equilibrium properties of 
this system are influenced by the transparency of the SN 
interfaces and by the presence of a spin-flip mechanism in 
the diffusive wire. In particular, we want to study (i) the 




FIG. 1: (Color online) Schematic representation of the sys- 
tem: a metallic diffusive wire (N) is connected at its ends to 
superconducting reservoirs Si and S r . Eventually, a supercur- 
rent / may circulate along the SNS junction. 



equilibrium density of states (DOS) in the normal wire 
and (ii) the supercurrent in the SNS system when a su- 
perconducting phase difference <fi is established between 
the electrodes. 

In order to describe these properties we use the qua- 
siclassical theory of superconductivity in the diffusive 
limit 3 ' where the mean free path is much smaller 
than the superconducting coherence length in the nor- 
mal metal, £ = y/ hD j A. This theory is formulated in 
terms of momentum averaged Green's function G(R, e), 
which depend on position R and an energy argument e, 
since we shall only deal with stationary situations. This 
propagator is a 4 x 4 matrix in Keldysh space (indicated 
by an inverted caret), where each entry is a 2 x 2 matrix 
in electron-hole space (indicated by a caret) 
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The general definitions of the different functions can 
be found in Ref. [34(. The Green's functions for the 
left (1) and right (r) leads can be written as Gj(e) — 

e -^ 3 ^/2SQ ( £ ) e *f 3 /2fi i where ^. ig the phase Q f the Qr _ 

der parameter of the electrode j = I, r. Here, Go(e) is the 
equilibrium bulk Green's function of a BCS superconduc- 
tor. Notice that, since we shall only consider equilibrium 
situations, the Keldysh component of G(R, e) can be ex- 
pressed in terms of the retarded and advanced compo- 
nents as G K = (G R - G A ) tanh(/3e/2), where (3 = l/k B T 
is the inverse of the temperature. 

The propagator G(R, e) satisfies the stationary Usadel 
equation, which in the N region reads 3 ^ 
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[f 3 Gf 3 ,G] +e[f 3 ,G] =0, (2) 



where fs is proportional to the unit matrix in Keldysh 
space and equal to the Pauli matrix T3 in electron-hole 
space. Equation |(5J) is supplemented by the normaliza- 



tion condition G~ 



-7T 1. In the previous equation, 



r s f is the scattering time associated to spin-flip (mag- 
netic) impurities or related pair-breaking mechanisms. 
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For instance, as it has been shown in Refs. ( lllJ 3 ll | . if 
the normal wire is a thin film and its width W does not 
exceed £, the effect of a perpendicular magnetic field H 
can be described with an effective spin-flip scattering rate 
r s/ = h/r sf = De 2 H 2 W 2 /(6h). 

In order to solve numerically the Usadcl equation it is 
convenient to use the so-called Riccati parametrization, 36 
which accounts automatically for the normalization con- 
dition. In this method and for spin-singlet supercon- 
ductors, the retarded and advanced Green's functions 
are parametrized in terms of two coherent functions 
7 
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(R, e) and 7 K ' A (R, e) as follows 
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with the "normalization matrices" 

R,A_ / (1 + 7^7^)-! 
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Some of these functions are related by fundamental 
symmetries (particlc-hole, retarded-advanced) like 

7 A (R, e) = - [7 fl (R, e)] * ; 7 A (R, e) = -J R (R, -*m 

Therefore, we just have to determine, for instance, the 
retarded functions. Using their definition in Eq. © and 
the Usadel equation ([2]), one can obtain the following 
transport equations for these functions in the normal wire 
region" 3 
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Here, x is the dimensionless coordinate which describes 
the position along the N wire and ranges from (left lead) 
to 1 (right lead). The expressions for T R , G R , T R and Q R 
are obtained by comparing Eq. ([T]) with Eq. (J3J. Notice 
that Eqs. (|5l6p couple the functions with and without 
tilde. This means in practice that, in general, one has to 
solve Eqs. (|5|6p simultaneously. 

Now, we have to provide the boundary conditions for 
the Eqs. (|5l6p . Let us first remind that for ideal interfaces 
(perfect transparency) such conditions at the ends of the 
N wire result from the continuity of the Green's functions 
over the SN interfaces: 

7%) = 7o%) ; 7%) = -7o%) 
^(e) = e-^ R (e-eV) ; j R (e) = -e^ 7( ?(e + eV%7) 

where 7; fl (e) = r y R (x — 0, e) and -f R (e) = 7^(2; = 1, e), 
and the same for the coherent function with tilde. Here, 



7o fl (e) = -A/{e R + ^ V /A 2 - (e«) 2 }, where e R = e + i0+. 
Finally, <p is the eventual phase difference between the 
two superconducting reservoirs, which we assume to be 
applied in the right electrode. 

For non-ideal interfaces one has to use the more gen- 
eral boundary conditions derived in Refs. [38l|3^ |. These 
conditions for an spin-conserving interface are expressed 
in terms of the Green's functions as follows 



G?d x & 3 = 
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Here, refers to the Keldysh-Green's function on 

side /3(a) of the interface, Go = 2e 2 /h is the quantum of 
conductance, Gn is the conductance of the normal wire 
and Ti are the different transmission coefficients charac- 
terizing the interface. In general, one would need the 
whole set {t^}, but since one does not have access to this 
information we adopt here a practical point of view. We 
assume that all the N interface open channels have the 
same transmission r and define Gb = GqNt as the con- 
ductance of the barrier. Thus, the two S-N interfaces 
will be characterized by two quantities, namely the bar- 
rier conductance Gb and the transmission r, and our 
starting point for the boundary conditions will be 



2tt 2 [G^,G a ] 



4vr 2 - T ({&,G a } +27T 2 )' 



(9) 



where we have defined the ratio r = Gn /Gb ■ In this lan- 
guage, an ideal interface is characterized by r = and a 
tunnel contact is described by t < 1. In what follows, 
unless the opposite is explicitly stated, we shall assume a 
symmetric situation with two identical interfaces. In the 
literature the so-called Kupriyanov-Lukichcv 23 boundary 
conditions are often used. These conditions can be ob- 
tained from Eq. ([9]) by removing the term proportional 
to t in the denominator. Such approximation is valid 
strictly speaking for the case of tunnel junctions (t « 1) 
and it turns out to be very good for highly transparent 
interfaces (r -C 1). 

The next step is to express these boundary conditions 
directly in terms of the coherent functions. Substituting 
the definitions of Eq. © into Eq. ^ and after straight- 
forward algebra, one obtains the following boundary con- 
ditions for the parameterizing functions 



(l-7^)7^-(l-7^)7^ 
(1 + 7#7|)(1 + 1 R 1 R ) ~ r(7 R ~ 7£)(7* - tf) 



U0) 



where the minus sign is for the left interface and the plus 
sign for the right one. The boundary conditions for 7^ 
can be obtained from Eq. (|10p exchanging the quantities 
without tilde by the corresponding ones with tilde and 
vice versa. These equations establish a relation between 
the functions and their derivatives evaluated on the side 
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of the interface inside the N wire (/?) and the correspond- 
ing functions evaluated on the side of the interface inside 
the reservoir (a), which are given by Eq. 10. 

In the limit of weak proximity effect, Eqs. (|5l6p can 
solved analytically, as we discuss in Appendix A. How- 
ever, in general they have to be solved numerically. These 
are second order differential equations with boundary 
conditions relating the functions and their derivatives 
in the two SN interfaces. This is a typical two point 
boundary value problem that we solve numerically using 
the so-called relaxation method as described in Ref. • 
We want to point out that the Riccati parametrization 
facilitates the numerical calculations because the coher- 
ent functions are smooth and bounded. Moreover, this 
parametrization is also well suited for time-dependent 
problems, as we have shown in Ref. [41(. 

To end this section we discuss the formula for the su- 
percurrcnt. The electrical current can expressed in terms 
of the Usadcl Green's functions as22 



87r 2 e 



r deTr{f 3 [G^G] K ( e )} 

J — OO 



(11) 



Combining this expression with fundamental symmetries 
of the Green's functions and using the fact that we only 
address equilibrium situations, we can write the super- 
current as 



Gn 

e 



deS(e) tanh(/3e/2), 



(12) 



where S = (l/4:TT 2 )Re{Tr(T 3 G R d x G R )} is the spectral 
supercurrent. 

III. LOCAL DENSITY OF STATES 

In equilibrium, the most basic quantity that reflects 
the proximity effect in the N wire is the local den- 
sity of states (DOS), which is defined as DOS(x,e) = 
— Im {Q R (x, e)} /it. This quantity can in principle be 
measured with a tunneling probe electrode as in Ref. [|| , 
or with a scanning tunneling microscope (STM) as in 
Ref. [4|. In this section we analyze the local DOS in the 
normal wire in different situations. 

Let us start discussing a situation where there is no 
phase difference between the superconducting reservoirs 
((f) = 0). In this case, one can show that the relation 
7 fl (e) = — "f R (e) holds. Thus, one only needs to solve 
Eq. §5§ for the coherent function 7 fl (e). In Fig.[2]we show 
an example of the local DOS in the middle of a normal 
wire (x = 0.5) of length L = 2£ without spin-flip scatter- 
ing (T s f — 0). The most prominent feature is the appear- 
ance of a minigap A g in the spectrum, which for perfect 
transparency scales with the Thoulcss energy roughly as 
A 9 ~ 3.1et in the long junction limit (L 3> £). Let us re- 
mind that the minigap is the same along the normal wire, 
although the exact DOS depends on the position. The 
existence of a minigap in a diffusive normal metal in con- 
tact with a superconductor was discussed by McMillan 9 




FIG. 2: (Color online) Density of states of a SNS junction as 
a function of energy in the middle of a wire of length L — 2£ 
without spin-flip scattering (F s f — 0). The SN interfaces 
are assumed to be identical and there is no phase difference 
between the S electrodes. In panel (a) the different curves 
correspond to different values of the ratio r = Gn/Gb and 
a transmission r = 1, while in panel (b) they correspond to 
different values of the transmission r for a ratio r — 1.0. 



within a tunneling model, where the normal region was 
a thin layer. In more recent years, the minigap has been 
extensively studied in various hybrid diffusive SN and 
SNS structures! 10 ' 11 ! 12 As one can see in Fig. HJa), this 
minigap diminishes progressively as the ratio r increases, 
i.e. as the interface becomes more opaque. For this par- 
ticular length, we find that the minigap decays with the 
interface parameter r as A g /A ~ 0.14/r for r > 1 (sec 
the fit to our numerical data in Fig. [T31 Appendix B). 

In Fig. HJb) we illustrate the effect of the transmission 
coefficient r in the local DOS for a ratio r = 1.0. Notice 
that the minigap is only slightly reduced as r decreases, 
while the features around A g become more pronounced. 
The effect of a transmission smaller than one is much 
more pronounced for larger values of r, i.e r 3> 1, while 
for values r < 1 it is rather insensitive to the value of r. 

In Fig. [3] we present a detailed study of the decay of 
the minigap as a function of the wire length for differ- 
ent values of the interface resistance and r = 1.0. We 
have normalized the minigap A g with the Thouless en- 
ergy to show explicitly that in the long wire limit A g sim- 
ply scales with this energy. In this limit (A/ey — > oo) we 
were able to fit accurately the decay of the minigap with 
the ratio r with the function A g /e T = 0.64/(0.20 + r) 
(see Fig. [T3]in Appendix B). In the opposite case of a 
short junction, i.e. when L <^ £, the minigap is of the 
order of A for perfect transparency, while it is given by 
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FIG. 3: (Color online) Minigap A 9 of a SNS junction as a 
function of the length of the N wire for different values of 
r — Gn/Gb and r s / = 0. The contact is assumed to be 
symmetric and the transmission is set to r = 1. Notice that 
is A g is normalized by the Thouless energy ex- 



A g - e T /2r in the limit of r > lM 

Let us now study how the density of states is mod- 
ified when there is a finite phase difference <f> between 
the leads, i.e. in the presence of a supercurrent. fn this 
discussion we shall assume that T s f = 0. Considering 
ideal interfaces, Zhou et al^ showed theoretically that 
the minigap decreases monotonically as the phase differ- 
ence increases and it closes completely when <p = it. In 
Fig. [4] we show two examples for L = 2£ of how the DOS 
in the middle of the wire evolves with the phase <fi for 
perfect transparency and r = 1.0. Notice that for finite 
r the qualitative behavior of the minigap is very similar. 
Indeed, a detailed study shows that, if the minigap is nor- 
malized by its value at <fi = 0, its phase dependence does 
not change significantly with the interface resistance. No- 
tice, however, that the features in the DOS around the 
minigap can be clearly different, as Fig. [5] exemplifies. 

Now we turn to the analysis of the influence of spin-flip 
scattering in the local DOS. Belzig et alJ^- showed that 
the minigap of an SN structure is reduced in the presence 
of a spin-flip mechanism and vanishes for large values of 
r s /. Different author a 24 ' 29 i 30 have studied the effect of 
magnetic impurities in the transport of SN structures and 
found that the Thouless energy is the scale that controls 
the effect of spin-flip on the proximity effect. In par- 
ticular, Crouzy et al}^ have shown analytically that in 
the long junction limit of an SNS structure, the minigap 
closes at a critical value of w 4.96et- 

Fig. displays the local DOS in the middle of a normal 
wire of length L = 10£ (er = 0.01A) for different values 
of the spin-flip rate T s f. The upper panel shows the case 
of ideal interfaces, while the lower one contains the results 
for a ratio r = 1.0. One can see how the minigap is pro- 
gressively reduced as F s f increases and finally vanishes. 
For r — (perfect interfaces) we find numerically that 
the gap closes at r^y w A.9ex in very good agreement 
with the long junction limit mentioned above^ For the 
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FIG. 4: (Color online) Density of states of a SNS junction as 
a function of energy in the middle of a wire of length L — 2£ 
(r s / = 0) for different values of the superconducting phase 
difference 4>. The interfaces are identical and characterized 
by a ratio r = Gn/Gb = 0.0 in panel (a) and r = 1.0 in panel 
(b). In both cases r = 1. 



case r — 1.0 this critical value is Tj^ w 0.8er- This means 
that it is reduced by approximately a factor 6, which is 
the same reduction factor obtained for the minigap (see 
Fig. [3]). This indicates that at finite transmission the rel- 
evant scale for the proximity effect is the minigap rather 
than the Thouless energy. This will become even clearer 
in the analysis of the supercurrent in the next section. 



IV. SUPERCURRENT 

As mentioned in the introduction, the supercurent in 
diffusive SNS junctions has been the subject of numer- 
ous theoretical and experimental studies. In particular, 
from the theory side, the results for the critical current 
for ideal interfaces and without spin-flip scattering are 
summarized in Refs. [6ll25j. The critical current in SNS 
junctions with partially transparent interfaces was dis- 
cussed in Ref. [231 ] using the boundary conditions devel- 
oped in the same reference. More recently, Heikkila et 
al£^ studied the reduction of the zero-temperature crit- 
ical current with the interface resistance considering a 
disordered interface. 

In this section we shall discuss how both the super- 
current and the critical current are modified by a finite 
transparency of the interfaces. To be precise, we shall 
investigate both the current-phase relationship and the 
temperature dependence of the critical current. More- 
over, we shall study in detail the effect of spin- flip scat- 
tering in the critical current, which to our knowledge has 
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FIG. 5: (Color online) (a) Density of states of a SNS junction 
as a function of energy in the middle of a wire of length L = 
10£ for r = and r = 1.0. The different curves correspond to 
different values of the spin-flip scattering rate r s /. (b) The 
same but for r = 1.0. Notice that both the energy and the 
r s / are in units of the Thouless energy. 




FIG. 6: (Color online) Zero-temperature supercurrent-phase 
relation of a diffusive SNS junction for L = 8£ (A = 64er) 
for different values of the ratio r = Gn/Gb and r = 1.0. In 
panel (a) we show the results for the product eRI/A, where 
R is the normal-state resistance, while in panel (b) we have 
normalized the different curves by the critical current Ic- 



not been discussed before in the literature. This analy- 
sis is very relevant from the experimental point of view, 
since it might describe the supercurrent in the presence 
of a magnetic field, as explain in the previous section. 
Finally, let us remind that the results of this section are 
complemented with Appendix A, where we study ana- 
lytically the supercurrent for the case of low transparent 
interfaces (r>l and t<1). 

We start our discussion by analyzing the current-phase 
relation in the absence of spin-flip scattering. In Fig. [5] 
we show this relationship at zero temperature for a wire 
of L = 8£ for different values of the ratio r. As it can 
be seen, the supercurrent is a non-sinusoidal function 
of the phase difference, which reaches its maximum at 
(\> k 1.277r/2, almost irrespectively of the value of r. For 
the ideal case (r = 0) this result agrees with the previ- 
ous results reported in the literature.— It is important to 
stress that in this figure and in what follows we normalize 
the current with the total resistance in the normal state, 
i?, which includes the contributions of both the diffusive 
wire and the interfaces. For a symmetric junction this 
resistance can be expressed in terms of the ratio r as 
R = (1 + 2r)/G N . 

Notice that, as one can see in Fig. UJb), when the su- 
percurrent is normalized by the critical current Ic, the 
different results almost collapse into a single curve. At 



a first glance, this result seems to suggest that the in- 
terface transparency just enters as a reduction prefactor 
in the expression of the critical current. However, as we 
discuss in the next paragraph, this is clearly not the case 
at finite temperature. 

Let us now turn to the analysis of the temperature de- 
pendence of the critical current Ic- In Fig. [7] we show 
this dependence for a wire of length L = 8£ and dif- 
ferent values of r. Notice that the temperature is nor- 
malized with the Thouless energy. The main conclusion 
that can be extracted from these results is that the crit- 
ical current decays faster with temperature as the inter- 
face resistance increases. Moreover, the saturation re- 
gion at low temperatures in which the critical current 
is almost constant shrinks as the interface resistance in- 
creases. For ideal interfaces (r = 0) this region corre- 
sponds, roughly speaking, to the range k^T < ct, while 
for finite r it corresponds to k^T < A g . This illus- 
trates the fact that is the minigap the scale that controls 
the magnitude of the supercurrent at arbitrary trans- 
parency. The faster decay for partially transparent inter- 
faces can be confirmed analytically in the limit of very 
long junctions (ey/A — * 0). In this case and for per- 
fectly transparent interfaces (r — 0), one finds a critical 
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FIG. 7: (Color online) (a) Critical current of a diffusive SNS 
junction as a function of the temperature for L = 8£ and dif- 
ferent values of the ratio r = Gn/Gb and r = 1.0. (b) The 
same but normalized by the critical current at zero tempera- 
ture. 
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4>/7T 

FIG. 8: (Color online) (a) Zero-temperature critical cur- 
rent of a diffusive SNS system with r = 1 as a function of 
the parameter r — Gn/Gb for different lengths of the nor- 
mal wire, (b) Current-phase relation at zero temperature 
for an asymmetric junction with L — 2£ and r = 1. The 
asymmetry parameter n is defined as k = 1 — tr/vl with 
f'R < ft and ft + fa = tlr = const. Here Tlr = 2 and 
ft = 1, 1.5, 1.7, 1.8, 1.9, 1.95, 1.99999 



current that decays as Ic oc (hsT/ex) 3 ^ exp(— L/Lt), 
where Lt — \JhD /2irkBT is the thermal length (see 
Refs. [6,22]]). In the opposite case of opaque interfaces 
(r ^> 1), the result of Appendix A indicates that the crit- 
ical current decays as Ic oc (hsT/ex) 1 ^ exp(— L/Lt)- 

The decay of the zero-temperature critical current with 
the interface resistance is examined systematically in 
Fig. [U^a) for different wire lengths and fixed transmis- 
sion t = 1. After normalizing the curves by the resis- 
tance in the normal state we find in the limit of very 
short wires (L « Q a maximal critical current at finite 
f before it slowly decays for large interface resistances. 
Thus, eRIc ^ A in the whole parameter space. For wires 
with L > £ we find a monotonic decay of eRIc/ A with 
increasing r. Then, for A/ct —> oo the energy scale of 
the critical current for large ratios r is determined by an 
effective Thouless energy £T,cff/eT ~ A r B /{C + r). For 
instance, when r > 10 we can fit the decay of the eRIc 
product for the special case of a wire with L = 9£ with 
the help oieRI c /e T = 5.13 f a29 /(0.22 + f) (see Fig. US). 
Here a fitting curve with B = would be proportional 
to the minigap but would only give a rough estimate of 
eRIc{r). So far we do not have a good explanation of 
the factor r B and the numerical value of B. 

The lower panel of Fig. [5] shows the current-phase re- 
lation at zero-temperature for a junction with L = 2£, 



t = 1 and asymmetric barriers as a function of the asym- 
metry parameter k = 1 — ft/r_R that fulfills r l + f _r = 
tlr = const. The critical current shows an enhancement 
for larger asymmetries while the phase difference moves 
towards 7r as k increases. By modeling the diffusive SNS 
junction as a point contact and averaging the current 
through the different channels over the bimodal distribu- 
tion for diffusive systems, one can understand this trend 
with the help of the Kirchhoff rules and the set of possi- 
ble shapes of the current-phase relation in this regime.— 
Furthermore, the formulas of the Appendix A can be gen- 
eralized to the asymmetric case. Then, the eRIc prod- 
uct, Eq. (|A9j) . is proportional to (ft + ffl)/(rt,r_R) what 
is in agreement with our numerical results. 

Let us now discuss the influence of a spin-flip mecha- 
nism in the supercurrent. As explained above, the spin- 
flip scattering may be due to paramagnetic impurities 
and in this case r s / is proportional to the impurity con- 
centration, or it may be caused by a magnetic field and 
in this case T s f is proportional to the square of the field. 
Indeed, this second possibility is much more interesting 
since it offers a natural way to control the strength of 
the spin-flip scattering and, in this sense, it is also more 
relevant from the experimental point of view. 
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FIG. 9: (Color online) (a) Zero-temperature critical current of 
a diffusive SNS as a function of the spin-flip rate for different 
lengths of the normal wire and ideal interfaces (r = 0). (b) 
The same but the critical current is now normalized by the 
zero spin-flip rate value and F s f is normalized by the Thouless 
energy of the wires. Notice that all the curves collapse into a 
single one that can be approximately described by the Gaus- 
sian function Ic/Ic(Xsf = 0) = exp(— 0.145r s / /er) shown as 
a back solid line. 



Fig. [9] displays the zero-temperature critical current as 
a function of the spin-flip rate T s f for different values of 
the wire length and ideal interfaces (r — 0). The reason 
for plotting the current as a function of the square root of 
the rate is that this plot can be seen as the magnetic field 
dependence of the critical current when the normal wire 
is a thin film. It is important to remark that in these cal- 
culations we assume that the order parameter in the leads 
is not affected by the spin-flip mechanism (such an effect 
can be trivially included) . As one can see in Fig. [9] the 
spin- flip mechanism causes a decay of the critical current. 
It is well-know n 27 ' 28 that the spin-flip scattering acts as a 
pair-breaking mechanism for the Cooper pairs that pen- 
etrate in the normal wire. Such scattering introduces 
a new relevant length scale in the problem, namely the 
spin-flip length L s f = y/hD/2T s f. When this length be- 
comes smaller than the length of the system and the ther- 
mal length, it dominates the decay of the supercurrent. 
As we show in Fig. HJb), when Iq its normalized by its 
value in the absence of spin-flip rate, its decay with T s f 
becomes universal for relatively long wires. Such decay 




FIG. 10: (Color online) (a) Zero-temperature critical current 
for a wire length L = 10£ as a function of the spin-flip rate 
for different values of the interface resistance and r = 0. (b) 
The same as in the upper panel, but the critical current is 
normalized by its value at F s f = 0. 



can be phenomenologically fitted with a Gaussian func- 
tion Ic/Ic(T s f = 0) = exp(— 0.145r s //ey), as demon- 
strated in Fig. [9jb) . The analysis detailed in Appendix 
A suggests that, in the low transparent regime, the decay 
follows a law of the type Iq oc (er/2r s j) 1 / 2 exp(— L/L s t) 
at finite temperature, which numerically is similar to the 
Gaussian function above. 

On the other hand, as one can see in Fig. [UJ there is still 
a non-negligible supercurrent even when minigap is com- 
pletely closed, i.e. when T s f > 5er- This phenomenon 
in the proximity structure considered here is the equiva- 
lent of the well-known gapless superconductivity in bulk 
superconductors ! 27 ' 28 

In order to understand the role of the interface trans- 
parency in the decay of the critical current as a function 
of r s / , we present in Fig. [10] the results for Iq for a wire 
of length L = 10£ for different values of the ratio r. As 
it can be seen in particular in Fig. llOf b). the critical cur- 
rent decays faster as the interface resistance increases. 
This fact illustrates again that the most relevant energy 
scale at finite transparency is the minigap rather than 
the Thouless energy. 

Finally, to complete the discussion of the role of the 
spin-flip scattering, let us now describe what happens at 
finite temperatures. In Fig. QT] one can see the critical 
current for a wire length L — 10£ as a function of the 
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T = 0.0 




FIG. 11: (Color online) (a) Critical current for a wire length 
L = 10£ as a function of the spin-flip rate for different values 
of the temperature (in unit of the critical temperature, To), 
(b) The same as in the upper panel, but the critical current is 
normalized by its value at T s f — 0. The inset show the tem- 
perature dependence of the constant a used in the Gaussian 
fits: I c /Ic(Tsf = 0) = exp(-ar s/ /e T ). 



rate T s f for different values of the temperature. The 
main conclusion is that by increasing the temperature 
the decay of the critical current becomes slower. Such a 
trend can be understood with the help of the result of 
Appendix A. 



V. CONCLUSIONS 

With the advances in the fabrication techniques of su- 
perconducting hybrid structures and the development of 
local measuring probes, it is now possible to explore the 
proximity effect in diffusive metallic nanostructures in 
great detail. In this sense, it is highly desirable from the 
theory side to elucidate the role of ingredients usually 
present in experiments such as partially transmissive in- 
terfaces and pair-breaking mechanisms. With this idea 
in mind, we have presented in this work a detailed anal- 
ysis of the density of states and supercurrent in diffusive 
SNS junctions. In particular, we have studied the influ- 
ence in these two equilibrium properties of an arbitrary 
transmission of the interfaces and spin-flip scattering in 
the normal wire. Our analysis is based on the quasiclassi- 



cal theory for diffusive superconductors (Usadel theory), 
supplemented by the boundary conditions put forward 
by Nazarov^ 

With respect to the local density of states, we have 
shown that the minigap that appears in the normal wire 
is very sensitive to the interface transmission both in the 
absence and in the presence of a supercurrent in the sys- 
tem. Moreover, we have shown that the minigap closes 
when the energy rate that describes the spin-flip scatter- 
ing is a few times larger than the minigap in the absence 
of this type of scattering. This fact nicely illustrates that 
the minigap is indeed the relevant energy scale for the 
proximity effect for non-ideal interfaces. 

Turning to the analysis of the supercurrent, we have 
shown that both the magnitude and temperature depen- 
dence of the critical current depend crucially on the inter- 
face resistance. In particular, the critical current decays 
faster with temperature as the interface resistance in- 
creases. Moreover, we have studied how the existence of 
spin-flip scattering in the normal wire diminishes the su- 
percurrent and identified the relevant energy and length 
scales for its decay. In particular, we have shown that 
a supercurrent can still flow when the minigap is com- 
pletely closed, which is the analogous in proximity struc- 
tures of the well-known gapless superconductivity in bulk 
samples i 27 i 28 This prediction can be tested experimen- 
tally by using an external magnetic field, as long as the 
width of the normal wire is smaller or comparable to the 
superconducting coherence length ad it is made of a thin 
film^i 
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APPENDIX A: LINEARIZED EQUATIONS 

In the limit of very low transparent interfaces (r ^> 1 
and r < 1), the supercurrent can be computed analyt- 
ically by linearizing the Usadel equations- 23 In this ap- 
pendix we describe how this can be done within the for- 
malism presented in section II. 

Assuming that the proximity effect in the normal wire 
is weak, the coherent functions are small and the retarded 
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and advanced Green's functions can be approximated by 
(see Eq. © 

G RA « ( ^ ) . (Al) 

Here, the coherent functions 7^ and j R fulfill the lin- 
earized version of Eqs. (|5l6p , which reduce to 



(A2) 
(A3) 



Notice that now the equations for ^ R and 7^ are uncou- 
pled and have an identical form. 

The boundary conditions for the previous equations 
are obtained by linearizing Eq. (|10[) in the following way 





7a 


T R 

J a 


1 - 


f 7a 7a 


2tu 






f R 


1 - 


f 7a 7a 


2wi 



(A4) 



where the minus sign is for the left interface and the 
plus sign for the right one. Here, J- R and F R are the 
anomalous Green's functions of the corresponding super- 
conducting lead a — l,r. 

The solution of Eq. (|A2[) with the boundary conditions 
of Eq. (|A4[) can be written as 



7 fl (x) = A R e iXx + B R e~ iXx , 



(A6) 



where A 2 = 2{ie — T s f)/eT and the constants A R and B R 
can be expressed as 



A 1 



B 1 



AnirX sin A 
1 

Airir A sin A 



(T R +T R e-^) 



{T R + T R ^ 



The solution for the function j R is obtained from the 
solution for j R by replacing the functions T R { by T R V 

After linearizing the expression of Eq. (fT2|) . the super- 
current can be written as 

I=^rf deRe {j R d x7 R - 7 R d x j R } tanh (^j . 

(A7) 

Using the solutions for -f R and j R , it is straightfor- 
ward to show that the supercurrent-phase relation can 
be written as 

/ = sin(0) ^ de Re tanh (%) , 

e-R z r z J-00 y2i\sm\ J \ 2 / 

(A8) 

where J- R is the bulk anomalous Green's function with- 
out including the superconducting phase. This integral 
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FIG. 12: (Color online) Comparison between the linearized 
result of Eq. (|A9|) and the exact numerical result for the tem- 
perature dependence of the critical current for a wire with 
£T = 75 A in the absence of spin- flip scattering (r s / = 0). 
The different curves correspond to different values of the ra- 
tio r = Gn/Gb and in the exact result we have used r = 0. 



can be done analytically and the result for the critical 
current is 



A 2 /(A 2 + c 2 ) 



where uj n = (2n + l)7rfcBT. Here, we have used R = 
(1 + 2r)/G^ ~ 2r/G*N- If in particular the temperature 
is just a few times larger than the Thouless energy, one 
just need to keep the first term (n — 0) in the previous 
expression. In the limit of an infinitely long wire this 
formula reduces to 



eRI c 



4nk B T 



exp(-L/L), 



(A10) 



where the effective length L = LtL s j / ^1 L\ + L s ^. Here, 

Lt = y/hD/2Trk-oT is the thermal length and L s f = 
y/HD/2T s f is the spin-flip length. 

In order to establish the range of validity of the ex- 
pression of Eq. (|A9[) we have compared this result with 
the full numerical solution of the non-linearized Usadel 
equations. An example of such a comparison for the tem- 
perature dependence of the critical current is presented 
in Fig.[l2]for a wire with ex = 75A. The different curves 
correspond to different values of the ratio r = Gn / Gb 
keeping always t = 0, which correspond to a tunnel junc- 
tion. Notice that the approximation of Eq. (|A9[) describes 
very well the exact results even for values of r very close 
to 1. 

In the limit of highly transparent interfaces, i.e. r <C 
1, the solution of Eqs. (|5l6p demands a more careful 
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treatment.? 2 - However, in the limit of high temperatures 
(k^T 3> A) one can obtain an analytical expression 
for the critical current using the linearized solution of 
Eq. (|A6[) and assuming that this function is continuous 
at the SN interfaces. This means in practice that the 
constants A R and B R appearing Eq. (|A6|) adopt now the 
following form 



1 



2i sm A 

B R = {rfe iX ~ 7^ 

2i sm A 



The rest of the calculation is identical and now the result 
for the critical current is 



eRI c = iTrkftT [ | <>xp( 



-L/L). 



(All) 



In the absence of spin-flip scattering this result re- 
produces the well-known result originally derived by 
Likharcv in Rcf. [21], which indicates that Ic <x 
(hsT/cr) 3 / 2 exp(— L/Lt)- It has been shown that this 
dependence describes a broad temperature rangep^ as 
long as ksT ^> €t- 
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FIG. 13: (Color online) Minigap and critical supercurrent as a 
function of r for the two lengths L = 2£ and L = 9£. The solid 
lines show the fitting-curves of the numerically calculated data 
(circles) in units of the Thouless energy for r ^> 1 . 



APPENDIX B: NUMERICAL FITS 

Fig. 1131 shows some of the numerical fits mentioned in 
sections III and IV. 
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